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Abstract 

We study the role of the tensor term of the Skyrme effective interactions on the spin-orbit split- 
tings in the N=82 isotones and Z=50 isotopes. The different role of the triplet-even and triplet-odd 
tensor forces is pointed out by analyzing the spin-orbit splittings in these nuclei. The experimental 
isospin dependence of these splittings cannot be described by Hartree-Fock calculations employing 
the usual Skyrme parametrizations, but is very well accounted for when the tensor interaction is 
introduced. The capability of the Skyrme forces to reproduce binding energies and charge radii in 
heavy nuclei is not destroyed by the introduction of the tensor term. Finally, we also discuss the 
effect of the tensor force on the centroid of the Gamow- Teller states. 
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Nuclei far from the stability valley open a new test ground for nuclear models. Recently, 
many experimental and theoretical efforts have been paid to study the structure and the 
reaction mechanisms in the nuclei near the drip lines. Modern radioactive nuclear beam 
facilities and experimental detector setups have revealed several unexpected effects in nuclei, 
such as the existence of haloes and skins [1] , the modifications of shell closures [2] and the 
pygmy dipole resonances [3, 4]. 

One of the current topics is the role of the tensor interactions on the shell evolution of 
nuclei far from the stability line. The role of the tensor interactions in the evolution of the 
single-particle states was first discussed within the Skyrme Hartree-Fock (HF) framework 
by Fl. Stancu et ai, almost thirty years ago [5]. However, serious attempts have never 
been devoted, until very recently, to the study of its effects on the evolution of the shell 
structure in heavy exotic nuclei. In fact, the Skyrme parameter sets which are widely used 
in nuclear structure calculations do not include the tensor contribution. This contribution 
was included only in the so-called Skyrme-Landau parametrizations of Ref. [6]. 

In the present paper, we discuss the isospin dependence of the shell structure (in partic- 
ular, the spin-orbit splitting) of the Z=50 isotopes and N=82 isotones. We use the HF plus 
Bardeen-Cooper-Schrieffer (BCS) approach, by employing a Skyrme parameter set plus the 
triplet-even and triplet-odd tensor zero-range tensor terms, which read 

v T = ^{[((7i • k')(a 2 ■ k') - ±((7! • a 2 )k >2 }5( ri - r 2 ) 
+ S(n - r 2 )[(<7i • k)(a 2 ■ k) - ^(<7i • a 2 )k 2 }} 

+ U{(a! ■ k')5(n - r 2 )(<7i • k) - l(<7i • a 2 )[k' ' 8(r 2 - r 2 )k]}. (1) 

In the above expression, the operator k = (Vi — V 2 )/2i acts on the right and k! = — (Vi — 
V 2 )/2i acts on the left. The coupling constants T and U denote the strength of the triplet- 
even and triplet-odd tensor interactions, respectively; we treat these coupling constants as 
free parameters in the following study. The tensor interactions (1) give contributions both 
to the binding energy and to the spin-orbit splitting, which are, respectively, quadratic and 
linear in the proton and neutron spin-orbit densities, 



jiiji + 1) - k{k + 1) - \ 



Rl(r). (2) 



In this expression q = 0(1) labels neutrons (protons), while where % = n,l,j runs over all 
states having the given q. The v 2 is the BCS occupation probability of each orbital and 
Ri(r) is the radial part of the wavefunction. It should be noticed that the exchange part 
of the central Skyrme interaction gives the same kind of contributions to the total energy 
density and spin-orbit splitting. The central exchange and tensor contributions to the energy 
density H are 

AH = \a{J 2 n + J 2 p ) + (3J n J p . (3) 
The spin-orbit potential is given by 
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where the first term on the r.h.s comes from the Skyrme spin-orbit interaction whereas 
the second term includes both the central exchange and the tensor contributions, that is, 
a = ac + olt and (3 = /3c + fir- The central exchange contributions are written in terms of 
the usual Skyrme parameters, 

ac = - h) ~ + t 2 x 2 ) 

o o 

(3c = -\{hxi + t 2 x 2 ). (5) 

o 

Basic definitions of all quantities derived from the Skyrme parameters can be found in 
Refs. [7, 8]. The tensor contributions are expressed as 

5 

a T = -U 

(3 T = ^(T + U). (6) 

The central exchange contributions (5) have been neglected when fitting most of the Skyrme 
parameter sets, and when performing most of the previous HF calculations. In this work, we 
employ the SLy5 parameter set [8] which has been fitted with the same protocol of the more 
widely used SLy4 set and should consequently have similar quality. In the case of SLy5, the 
central exchange contributions are included in the fit and we take them into account here. 

Except for the double-magic systems, we perform HF-BCS in order to take into account 
the pairing correlations. Our pairing force is a zero-range, density-dependent one, namely 

V = V (l- [^^j ) ' ~ *y. (7) 

The parameters of this force (that is, Vo=680 MeV-fm 3 , 7=1 and £>o=0.16 fm~ 3 ) have been 
fixed along the Z=50 isotopic chain in connection with the SLy4 set in Ref. [9]. Therefore, 
we employ the same force here both for the neutron and the proton pairing interactions 
in the 50-82 shell, neglecting small readjustments which could be made to account for the 
Coulomb anti-pairing effect in the case of protons. 

Before coming to a detailed analysis of our results, let us mention the important general 
features associated with the tensor and the central exchange contributions to the spin-orbit 
splitting. The first point concerns the A-dependence (or isospin dependence) of the the 
first and second terms in the r.h.s. of Eq. (4). Since the Skyrme spin-orbit contributuion 
(proportional to Wq) gives a value of the spin-orbit splitting which is linear in the derivatives 
of the proton and neutron densities, the associated mass number and isospin dependence is 
very moderate in heavy nuclei. On the other hand, the second term in Eq. (4) depends on 
the spin density J q which has a more peculiar behavior. J q gives essentially no contribution 
in the spin-saturated cases, but it increases linearly with the number of particles if only 
one of the spin-orbit partners is filled. The sign of the J q will change depending upon the 
quantum numbers of the orbitals which are progressively filled: that is, the orbital with 
j > = 1 + 1/2 gives a positive contribution to J q while the orbital with j < — I — 1/2 gives a 
negative contribution to J q . This must be kept in mind to understand the results which are 
discussed below. 

According to Ref. [5], the optimal parameters ax and (3t should be found in a triangle 
in the two dimensional Pt) plane, lying in the quadrant of negative ar and positive 
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Pt- At that time, the force SIII was used. As already mentioned, we wish to use here the 
Lyon forces which have been fitted using a more complete protocol (including the neutron 
matter equation of state) and have some better features like a more realistic value of the 
incompressibilty K^. Therefore, we have refitted the values of (ay, Pt) using the recent 
experimental data [10] for the single-particle states in the N=82 isotones and the Z=50 
isotopes. We have not tried to refit all the Skyrme parameters after including the tensor 
terms. This can be left for future work. However, we have checked that the binding energies 
and the r.m.s. charge radii of 132 Sn ( 208 Pb) change, respectively, by +0.65% and -0.17% 
(+0.46% and -0.11%) when we include the tensor force. The parameters we have chosen are 
(Xt = -170 MeV-fm 5 and fir = 100 MeV-fm 5 . We should mention that for the force SLy5, 
ac = 80.2 MeV-fm 5 and Pc — -48.9 MeV-fm 5 . The fact that we need significantly larger 
values of (ax, Pt) as compared to Ref. [5] can be adscribed to the fact that the effect of the 
central exchange terms, with (a c , Pc) having opposite sign as the (a-r, Pt) required by our 
fit, must be counterbalanced. 

In Fig. 1, the energy differences of the proton single-particle states, e(hu/2) — £(97/2), 
along the Z=50 isotopes are shown as a function of the neutron excess N-Z. The original 
SLy5 interaction fails to reproduce the experimental trend both qualitatively and quanti- 
tatively. Firstly, the energy differences obtained within HF-BCS are much larger than the 
empirical data. Secondly, starting from the double-magic 132 Sn isotope, the experimental 
data markedly decrease as the neutron excess decreases and reach about 0.5 MeV at the 
minimum value in 112 Sn. On the other hand, using the HF-BCS approach with SLy5 the 
result is qualitatively the opposite: the energy differences slightly increase as the neutron 
excess decreases, and there is a maximum around 120 Sn. We have studied also several other 
Skyrme parameter sets and found almost the same trends as that of SLy5. 

In the results displayed in Fig. 1, we can see a substantial improvement due to the 
introduction of the tensor interaction with (a T , /3t) = (-170,100) MeV-fm 5 . This choice gives 
a very nice agreement with the experimental data in the range 20 < (N-Z) < 32, both 
quantitatively and qualitatively. The HF+tensor results can be qualitatively understood 
by simple arguments, looking at Eq. (4). In the Z=50 core, only the proton g 9 / 2 orbital 
dominates the proton spin density J p (cf. Eq. (2)); consequently, with a negative value 
of a T) the spin-orbit potential (4) is enlarged in absolute value (notice that W is positive 
and the radial derivatives of the densities are negative), the values of the proton spin-orbit 
splittings are increased, and the energy difference e(h\i/2) — £{97/2) is reduced with respect 
to HF-BCS without tensor. This reduction is seen better around N-Z=20: in fact, 120 Sn 
is, to a good extent, spin-saturated as far as the neutrons are concerned so that one gets 
no contribution from J n . It should be noticed that the term in a does not give any isospin 
dependence to the spin-orbit potential for a fixed proton number, but only the term in P 
can be responsible for the isospin dependence. In a pure HF description, from N-Z=6 to 
14, the 517/2 neutron orbit is gradually filled and J n is reduced. Then, the positive value of 
Pt enlarges in absolute value the spin-orbit potential and increases the spin-orbit splitting, 
so that the energy difference e(hn/2) — £(97/2) becomes smaller. Because of pairing, this 
decrease is not so pronounced in the results of Fig. 1. Moreover, from N-Z=14 to 20, the 
s\/2 and ds/2 neutron orbits are occupied and in this region the spin density is not so much 
changed since the s\/2 orbital does not provide any contribution. Instead, for N-Z=20 to 
32, the hu/2 orbital is gradually filled. This gives a positive contribution to the spin-orbit 
potential (4) and the the spin-orbit splitting becomes smaller. e(hu/2) — £(97/2) consequently 
increases, and this effect is well pronounced in our theoretical results. The magnitude of 
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(3 determines the slope of the isospin dependence, so that a larger (3 would give a steeper 
slope. 

In Fig. 2, the energy difference £(^13/2)— ^(^9/2) for neutrons ouside the closed N=82 core is 
plotted as a function of the neutron excess. The notation is the same of the previous figure. 
Essentially, the same arguments already made in the previous paragraph can be applied 
in order to understand the results; simply, we should remind that the proton number is 
increasing from the right (where the last nucleus displayed is 132 Sn) to the left (where the 
first isotope plotted is 150 Er). The lg 7 / 2 and 2d 5 / 2 orbitals are rather close in energy, above 
the last occupied proton state lgg/2 of the Z=50 core, and their occupations are affected by 
the pairing correlations introduced by the BCS approximation. These two proton orbitals 
have opposite effect on the spin orbit potential (4). Because of its larger value of j, the 1^7/2 
orbital turns out to play a more important role on the spin-orbit potential, when the tensor 
interaction is included, in the nuclei with N-Z decreasing from 32 to 24. Accordingly, with 
positive (3t the neutron spin-orbit potential is enlarged in absolute value: the spin-orbit 
splitting is made larger for these isotones, so that the i 13 / 2 orbital is pushed down and the 
h 9 / 2 is pushed up. These changes make the energy gap £(^13/2) — £(^9/2) smaller for the 
nuclei from N-Z=32 ( 132 Sn) to N-Z=24 ( 140 Ce). Then, the occupation of the 2d 5/2 orbital 
reverses the trend around N-Z=22 ( 142 Nd). The theoretical trend remains the same until 
N-Z=14, since the effect of the 2g? 3 / 2 occupation is counterbalanced by the occupation of the 
l^-n/2 which is not much higher and enters the active BCS space. 

The role of the tensor interaction due to the /3t term is expected from the discussion made 
by J.M. Blatt and V.F. Weisskopf [11] for the deuteron. In Ref. [12] the same argument 
was also presented. The role of ay has not yet been examined in a quantitative way within 
the mean field calculations, as this term comes from the triplet-odd tensor interaction. The 
assessment of its role is new, since the triplet-odd tensor interaction was not included in 
the analysis of Refs. [11, 12]. Recently, Brown et al. [13] studied the tensor interactions 
in 132 Sn and 114 Sn, based on the parameter set Skx. They considered both positive and 
negative values of cxt in the HF calculations and they concluded that cxt < gives a better 
agreement with the experimental data. This result is consistent with the present systematic 
study of the single-particle states in the Z=50 isotopes and N=82 isotones, performed within 
the HF+BCS model. 

The effect of the tensor interactions can be tested on other single-particle states which are 
empirically known. In this work, we have also considered the relative position of the 2d%/2 
and 1 /in/2 neutron states in 132 Sn and 100 Sn. In the former case ( 132 Sn), experimentally 
the two states are the last occupied, with the 2d 3 /2 being less bound than lhn/2 by about 
240 keV (see Fig. 8 of [14]). Theoretically, all the mean field calculations with Skyrme or 
Gogny forces, as well as the relativistic mean field (RMF) ones, result with the opposite 
ordering (see Fig. 7 of [15]). In particular, with the SLy5 force employed here, the lh u / 2 
orbital is less bound by 1.76 MeV. The contribution of the added tensor force reduces this 
value to 0.67 MeV. It has to be noted that the the right position of the 2d 3 / 2 level may be 
rather important for the proper description of the low-lying dipole strength in 132 Sn. In fact, 
the results obtained using a Skyrme force in Ref. [16] show that this strength has basically 
single-particle character; but even in the calculation of Ref. [17], in which the low-lying dipole 
strength turns out to present a certain degree of collectivity, the energy position of the levels 
we have mentioned is relevant since the configurations involving the 2rf 3 / 2 hole contribute to 
about 50% of the low-lying collective state. In the nucleus 100 Sn, experimentally the 2d^/2 
level is more bound by 0.9 MeV (see Fig. 7 of [14]). In our HF calculation with the SLy5 
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force, the two levels have the right order but their energy difference is 2.13 MeV. It is quite 
satisfactory that this difference becomes 1.33 MeV after including the tensor interactions. 

We have already mentioned that total binding energies and charge radii of 132 Sn and 
208 Pb are not extremely sensitive to the tensor interactions. We have also checked in our 
work the effect of the tensor terms on the isotope shift of the Pb isotopes. Actually, this 
effect is small (and does not even have the correct sign to reproduce experiment). Thus, 
the tensor interactions is unable to produce the well-known empirical kink in the trend of 
the charge radii beyond 208 Pb, which instead results from the introduction of generalized 
spin-orbit functionals [18]. 

Single-particle energies, and other ground-state properties, are not the only observables 
which are affected by the inclusion of tensor interactions. There exist excited states which 
reflect very much the behavior of the spin-orbit splittings. One of them is the well-known 
Gamow- Teller resonance (GTR). We have made a simple estimate of the effect of the ten- 
sor interactions on the GTR centroid by using the sum rules. In charge-exchange RPA 
calculations, the following sum rules are satisfied [19], 

m_(0)-m + (0) = (0|[Ft,F]|0>, 

m_(l)+ro + (l) = (0|[Ft,[#,F]]|0>. (8) 

In the above expressions, m±(k) denotes the fc-th moment of the strength in the AT Z = ±1 
channel: in particular, we are considering the non energy-weighted sum rule m(0) and the 
energy-weighted sum rule m(l). The associated operators F and F^ act, respectively, in the 
AT Z = —1 and AT 2 = +1 channels. In the Gamow- Teller case, they read 

f = *(<)*-(<), 

i 

Ft = (9) 

i 

Moreover, in Eq. (8), H is the total Skyrme Hamiltonian and |0) is the HF ground state. 
In nuclei with neutron excess, the contributions associated with the AT Z = +1 channel are 
negligible and we can approximate the GT centroid as 

m_(0) (0|[Ft,F]| ) ■ 1 I 

The contribution from the tensor interaction to the GT centroid is obtained by replacing 
the total Hamiltonian H in the previous formula, with the two-body force of Eq. (1). The 
calculation of the ground-state expectation value of the double commutator [F\ [vt, F]} has 
been worked out and it gives the contribution of the tensor force to the GT centroid, 

An f 24 

AEgT = 9(N-Z) J drr \^ ~ 5a ^ J " - 12a ( J P + ( U ) 

We have evaluated this latter expression, for the Sn isotopes having N-Z larger than 20, 
by using our optimal (cxt, Pt) values of (-170, 100). The results are reported in Table I. 
The numbers are not small, but this should not be surprising since the shifts of the single- 
particle states displayed in the Figures can also be of the order of 1-2 MeV. In fact, the 
positive energy shift can be expected in 208 Pb as well and can be understood as follows. 
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The spin-orbit densities (2) receive contribution only from the i 13 / 2 orbital (in the case 
of neutrons) and from the hu/2 orbital (in the case of protons). From Eq. (4), one sees 
that, since ax is negative and \olt\ > \Pt\, the net effect is an increase of the spin-orbit 
splitting. Consequently, the excitation energy of the dominant unperturbed Gamow- Teller 
configuration, that is, 1^13/2 — > 7^11/2, is shifted upwards by the tensor correlations. With 
similar arguments we can understand the numbers reported in Table I, as we did for the 
values plotted in the previous Figures. Actually, we should remind that the analysis in 
terms of the sum rules is not able to tell whether the peak energy is affected as much as the 
centroid m(l)/m(0) since the main peak does not exhaust, as a rule, the whole strength. 

Accurate QRPA calculations of the Gamow- Teller and spin-dipole resonances are reported 
in Ref. [20]. The behavior of different Skyrme parameter sets, without the tensor contribu- 
tion, is critically discussed. In fact, no RPA or QRPA calculations including the two-body 
tensor force are presently available, for any kind of vibrational mode; accordingly, results 
obtained without the tensor interactions should be still kept as reference until a global refit 
of the Skyrme plus tensor parameters is carefully accomplished. The tensor force is not only 
expected to produce effect on the Gamow- Teller states. Other vibrational states (like the 
low-lying 2 + which in many systems, once more, reflects the spin-orbit splitting [21]) will be 
certainly affected. 

A further question for future work is the role of correlations beyond mean field. As 
discussed at length in Ref. [22], the coupling of single-particle states to vibrational states 
has the net effect of increasing the level density around the Fermi surface by about 30%, by 
shifting occupied and unoccupied states in opposite directions. Smaller effects are expected 
for the energy differences we are considering here since these differences involve pairs of 
states which are either occupied or unoccupied. The net shift may be of the order of few 
hundreds of keV as estimated from 132 Sn [23]. 

In conclusion, the present work has shed light on the necessity to include the tensor 
component in the Skyrme framework. The first attempts in this direction were focusing 
on the effect of the tensor force in magic nuclei but, as we have stressed, if the nuclei 
are spin-saturated the spin-orbit splittings are not affected at all by the tensor force. The 
experimental mesurement of the isospin dependence of single-particle energies has opened 
the possibility to fit the parameters of the zero-range effective tensor force we are employing. 
Our results show that the introduction of the tensor force can fairly well explain the isospin 
dependence of energy differences between single-particle proton states outside the Z=50 
core, and neutron states outside the N=82 core. We have not attempted to refit a Skyrme 
force by including the tensor contribution, but we have discussed, by using the case of the 
Gamow- Teller centroids, that excited state properties will also be affected by the tensor. An 
ambitious refitting program of Skyrme forces should therefore be undertaken and deformed 
systems should be considered as well [24]. This is left as a future prospect, together with 
the role of particle- vibration coupling in this context. 

Acknowledgments 

We would like to thank D. M. Brink for stimulating and enlightening discussions. One of 
us (G.C.) gratefully acknowledges the hospitality of the University of Aizu, where the present 
work has started. The work is supported in part by the Japanese Ministry of Education, 
Culture, Sports, Science and Technology by Grant-in- Aid for Scientific Research under the 



7 



program number (C (2)) 16540259. 



I. Tanihata et al, Phys. Rev. Lett. 55, 2676 (1985). 

P. G. Hansen and B. Jonson, Europhys. Lett. 4, 409 (1987). 

I. Tanihata, J. Phys. G22, 157 (1996). 

A. Ozawa, T. Kobayashi, T. Suzuki, K. Yoshida and I. Tanihata, Phys. Rev. Lett. 84, 5493 
(2000). 

A. Leistenschneider et al, Phys. Rev. Lett. 86, 5442 (2001). 
P. Adrich et al, Phys. Rev. Lett. 95, 132501 (2005). 

Fl. Stancu, D. M. Brink and H. Flocard, Phys. Lett. 68B, 108 (1977). 
K. F. Liu et al, Nucl. Phys. A534, 1, 25 and 48 (1991). 
T. H. R. Skyrme, Nucl. Phys. 9, 615 (1959). 

D. Vautherin and D. M. Brink, Phys. Rev. C5, 626 (1972). 

E. Chabanat, P. Bonche, P. Haensel, J. Meyer, R. Schaeffer, Nucl. Phys. A635, 231 (1998). 
S. Fracasso and G. Colo, Phys. Rev. C72, 064310 (2005). 

J. P. Schiffer et al, Phys. Rev. Lett. 92, 162501 (2004). 

J. M. Blatt and V. F. Weisskopf, Theoretical Nuclear Physics (Wiley, New York, 1952). 
T. Otsuka et al, Phys. Rev. Lett. 95, 232502 (2005). 

B. A. Brown et al, Phys. Rev. C, in press (2006). 

H. Grawe, M. Lewitowicz, Nucl. Phys. A693, 116 (2001). 

M. Bender, P.-H. Heenen, P.-G. Reinhard, Rev. Mod. Phys. 75, 121 (2003). 

D. Sarchi, P. F. Bortignon, G. Colo, Phys. Lett. 601, 27 (2004). 

D. Vretenar et al, Nucl. Phys. A692, 469 (2001). 

P.-G. Reinhard, H. Flocard, Nucl. Phys. 584, 467 (1995). 

N. Auerbach and A. Klein, Nucl. Phys. A395, 77 (1983). 

S. Fracasso and G. Colo, to be published. 

S. Peru, J. F. Berger, P. F. Bortignon, Eur. Phys. J. A26, 25 (2005). 

C. Mahaux, P. F. Bortignon, R. A. Broglia, C. H. Dasso, Phys. Rep. 120, 1 (1985). 
G. Colo and P. F. Bortignon, unpublished. 

T. Duguet, private communication. 



8 



TABLE I: The effect of the tensor force on the GT centroid, evaluated by means of Eq. (11) using 
the SLy5 parameter set and the same parameters of the tensor force which have been fitted on the 
experimental results for the single-particle states, namely (ay, /?r)=(-170,100) MeV-fm 5 . See the 
text for a discussion of these results. 



Nucleus 


AE GT [MeV] 


120 Sn 


1.49 


122 Sn 


1.55 


124 Sn 


1.74 


126 Sn 


1.99 


128 Sn 


2.21 


130 Sn 


2.48 


132 Sn 


2.64 




FIG. 1: Energy differences between the lh n / 2 and lg7/2 single-proton states along the Z=50 
isotopes. The calculations are performed without (crosses) and with (circles) the tensor term in 
the spin-orbit potential (4), on top of SLy5 (which includes the central exchange, or J 2 , terms). 
The experimental data are taken from ref. [10]. See the text for details. 



9 



Neutrons on N=82 core 
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FIG. 2: Energy differences between the li 13 / 2 and lhg/ 2 single-neutron states along the N=82 
isotones. The calculations are performed without (crosses) and with (circles) the tensor term in 
the spin-orbit potential (4), on top of SLy5 (which includes the central exchange, or J 2 , terms). 
The experimental data are taken from ref. [10]. See the text for details. 
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